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SUMMARY 


It is imperative that the vibration environment of future vehicles be predicted 
prior to design and development so that satisfactory design and test procedures can be 
established. These criteria are essential to the establishment of high reliability 
standards necessary for man rated vehicles. The methods and techniques presented 
herein allow adequate predictions of future vehicle environments. 

These techniques rely upon typical structural configurations which have been 
sufficiently defined by measured data. Subsequent statistical analyses describe the dy- 
namic characteristics of the structure with statistical certainty. Thus, with only a 
knowledge of the structural geometry and mass characteristics, the anticipated dynamic 
environment may be established. These techniques are applicable to all rocket vehicle 
structure including corrugated and sandwich skin construction. 

The predicted environments represent a statistical estimation since the reference 
spectra are established by statistical techniques. Consequently, the probability of the 
actual environment not exceeding the predicted environment of a future vehicle is 
established with a 97. 5 per cent confidence. This does not infer that the predicted 
environment will accurately correspond to a single measured environment. Certainly, 
some of the measured responses of a new vehicle will be significantly lower than the 
predicted. This is to be expected since the criterion is such that the prediction will 
envelope 97. 5 per cent of the situations. However, this problem is deviated somewhat 
by the techniques utilized of separating rocket vehicle structure into eight (8) basic 
categories - each possessing essentially similar dynamic characteristics. This reduces 
the variance about the mean so that the mode value (most likely to occur) is not greatly ' 
less than the higher confidence limits. Therefore, the 97. 5 per cent criterion may be 
used without the concern of over conservatism in regard to a specific problem. 


INTRODUCTION 


The problem 6f dynamic analysis has been greatly complicated in recent years by 
the nature of the exciting forces and subsequent effects. Dynamic environments caused 
by rocket vehicle operation may result in failures due to over-stress and/or fatigue. 
Consequently, dynamic environments must be specified prior to design and development 
so that adequate design and testing criteria can be established. The precise prediction 
of these environments is a highly complex criterion which may not be reached in the 
near future. However, based on measured data and a few simplifying assumptions, 
adequate predictions may be obtained provided the necessary assumption and limitations 
are realized by the user. 

Vibration originates primarily from four sources of excitation. These are: 

1. Mechanically induced vibration from rocket engine fluctuation which is trans- 
mitted throughout the vehicle structure. 


2. Acoustic pressures generated by rocket engine operation. 

3. Aerodynamic pressure created by boundary layer fluctuations. 

4. Self-generating machinery, etc. 

Restricting notation to vibrational power quantities, the total vibration at any 
point on the vehicle may be expressed as: 


P^ = P/ , + 
T mech 


P/ + P/ + P/ 


ac 


ae 


mach 


Where / denotes the vibrational power caused by the indicated source. The powers add 
arithmetically since vibrational power is proportional to the mean square cyclic re- 
sponse. An exact analysis of structural response would necessitate an accurate descrip- 
tion of each individual source and the manner in which they combine. In most cases only 
one source is the primary forcing function. Hence, the remaining sources may be con- 
sidered negligible in regard to the total dynamic response at any instant of time. How- 
ever, a technique is developed on page 12 which sufficiently solves the problem of com- 
bined driving functions. 

The current state of the art provides several methods for dynamic predictions. 
References 1 through 4 provide adequate techniques for estimating the vibrational re- 
sponse caused by various driving forces. There is widespread agreement that oscilla- 
tory forces are significant and should be accounted for in design and environmental test- 
ing. Conversely there is widespread disagreement as to how the dynamic environments 
should be established, specified, and utilized. A study of the wide differences in tech- 
niques and the amount of time required for analysis indicates the obvious necessity for 
standardization. The problem then is to present an adequate method, free of complex 


2 


calculations, for determining dynamic environments which may be used for design and 
test criteria. It is realized that the average engineer does not have available time to 
thoroughly understand and analyze the dynamic characteristics of each individual prob- 
lem. Consequently, a generalized method must be presented which sufficiently solves 
the average problem with a minimum amount of calculations. Also, it should be realized 
that each problem is unique and requires at least some minor calculations. Hence, an 
optimum point must be reached which meets the requirements of both simplicity and 
adequacy. The method presented in the subsequent section is currently considered to 
satisfy these requirements. The method allows prediction of vibratory environments 
with only a few simple computations and excludes all complexities which might cloud the 
issue. The prediction methods and techniques advocated herein have proved successful. 
Research programs are currently being conducted in order to more precisely define the 
degree of accuracy obtained from these methods and direct studies for future refinements. 
Based upon additional information regarding the dynamic response phenomena and labora- 
tory test results, the existing techniques may be slightly modified. However, the basic 
foundation upon which the methods were constructed will remain unchanged. 
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DESCRIPTION 


Rocket vehicle structure may be separated into two dynamic categories. These 

are: 


I. Structure susceptible to acoustic and aerodynamic pressures (i. e. , 

skin panels ) . 

II. Structure not susceptible to acoustic and aerodynamic pressures (i. e. , 
structural beams). 

The dynamic analysis of these two groups are obviously handled by two different tech- 
niques. Each will be explained in detail and supplemented by example problems. 

Type I 

First, consider Type I. This group of structure may be subdivided 

into three sections. 


1(a) — Skin panels. A panel is defined as a section of skin bounded 
by radial and longitudinal stiffeners. This type of structure is directly excited by im- 
pinging acoustic pressures. The direction normal to the panel face exhibits the most 
severe vibratory response; consequently, this is the direction under consideration when 
referring to panel vibrations. 

1(b) — Skin stiffeners (such as ring frames and stringers). This type 
of structure is not directly excited by acoustic forces but is driven by the motion of 
adjacent panels. Thus, the stiffeners are indirectly forced by impinging acoustic pres- 
sures. 


1(c) — Bulkheads. These skin segments form the upper and lower 
extremities of vehicle propellant tanks. The bulkheads are further subdivided into: 

1(c) -1: Forward bulkheads 

1(c) -2: Aft bulkheads 

The added mass of liquid loading greatly reduces the vibration amplitudes experienced 
by the aft bulkheads; consequently, the two bulkheads are treated separately. 
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Although. aerodynamic pressure has not been mentioned, it should be realized 
that structures susceptible to acoustic pressures are indeed sensitive to aerodynamic 
pressures. The nature of the exciting forces are essentially the same, resulting in a 
similar vibratory response. The maximum acoustic pressures generally occur during 
captive firing or the liftoff period of the vehicle. The maximum aerodynamic pressure 
occurs at some period during flight, generally subsequent to Mach one. 

The equation for predicting the vibration environment of acoustically (or aero- 
dynamically) susceptible structure is: 



where G = the vibration response of the new vehicle structure at a particular station 
number. The term G is the acceleration due to cyclic motion divided by 
the acceleration of gravity g 0 . Since rocket vehicle vibrations contain 
many frequencies the response magnitude (G) is specified in spectral form. 

G = the known vibration response of a reference vehicle structure. This value 
has been determined by measurements and is also presented on a spectral 
basis. 

t = the thickness of the skin associated with G 

r r 

= the skin weight density of the reference structure 

p = the impinging acoustic (or aerodynamic) pressure which is driving the 
reference structure 

t = the skin thickness 1 associated with G 

n n 

P n = the skin weight density associated with G n 

p = the acoustic (or aerodynamic) pressures impinging upon the new vehicle 
structure. This pressure must be predicted. 

F = a factor which accounts for the attenuation effects produced by incorporating 
additional mass into the existing system. 


1 For corrugated or sandwich structure this parameter is an equivalent thickness. The 
technique for determining this thickness is shown in Appendix I. 
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The above equation applies to random rms composite values or to sinusoidal 
values. For the mean square spectral value the total expression simply raised to the 
second power. 

Equation (1) is applicable to localized vibratory environments and is valid for 
all materials. It is invalid when considering large sections of vehicle structure (i. e. , 
entire cylindrical tank). However, the static loading of these large sections is the 
critical design factor and localized dynamics thereby produce only negligible effects. 
Further, the added mass effects of fluid loading are not considered for skin panels and 
stiffeners of the propellant tankage since the most severe environmental condition is the 
primary objective. It is assumed that all vehicle stages are captive-fired prior to 
launch. Thus, there will be a period of time in which the fluid level is below the point 
of interest resulting in extreme dynamic responses. To thoroughly understand the equa- 
tion, a derivation is provided in Appendix I. 

Referring again to equation (1) p n represents the maximum pressure impinging 
upon the vehicle at any time. Three conditions must be considered, one of which will 
result as the maximum. These are: 

1. Captive-firing' environments (applicable to boosters and upper stages) 

2. On-pad acoustic environments 

3. The period of maximum aerodynamic pressures (occurs subsequent to Mach 
one; therefore, the coupling of acoustic-aerodynamic pressures do not have to be 
considered) . 

A flow chart indicating the train of events necessary to determine the vibration 
environment of a new vehicle is shown in Figure 1. Sound pressure curves corres- 
ponding to all standard liquid rocket engines are provided in Figures 2 through 8. The 
acoustic environments produced as a result of rocket engine noise for the different types 
of engines presented herein are based on the assumption that the engine rated parameters 
of thrust, flow rate and exhaust velocity will be generated at sea level conditions. The 
acoustic environment is given in terms of overall rms (5 cps to 2000 cps) dynamic 
pressure, in pounds per square inch, and is plotted as a function of vehicle length. The 
vehicle length is measured in inches above the engine nozzle exit plane. The curves 
represent a single engine operation of each respective type of engine (i. e. , one M-l 
engine, one F-l engine, one J-2 engine, etc. ). To determine the overall acoustic pres- 
sure, at a given vehicle length, resulting from operation of two or more engines of any 
one particular type, the following equation can be used: 



where: N = number of engines of one given type 
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pj = the pressure produced, at a given vehicle length, from the operation of 
a single engine. (Figures 2 through 8) 



The above equation cannot be used to determine the resultant pressure from operation 
of two or more different types of engines, (i. e. , four M-l engines and six F-l engines 
on the same stage). 

The effect of rocket engine clustering is not well defined at this time; therefore, 
it was not explicitly treated. It is felt however, that the estimates of overall SPL ob- 
tained from Figures 2 through 8 are valid if the number of engines clustered are not 
great. It has been observed from measured acoustic data from clusters of one, two, 
four, and eight H-l engines which propel the booster of the Saturn I vehicle, that cluster- 
ing does not produce appreciable effects in peak frequency of the acoustic spectra. From 
model tests involving as many as 24 engines per cluster, only slight effects if any 
were exhibited in octave band frequency spectra. This model test should not be consider- 
ed conclusive. The measured data should, however, be considered as indicative since 
they do represent test results. Therefore, it is suggested that caution be taken for use 
of this technique for engine clusters of more than 16 engines. 

The psue do- acoustic pressures on the vehicle surface induced by the aerodynamic 
boundary layer during flight may, for some locations, exceed the engine sound pressures 
even during the on-pad condition. Figure 9 shows the maximum in-flight rms overall 
sound pressures anticipated for future space vehicles utilizing liquid propellant boostei 
stages. The words, "rms-overall" denote the total rms aerodynamic pressure in a 
spectral window from 5 cps to 2000 cps. Acoustic measurements from Saturn flights 
have shown a ratio of the overall aerodynamic oscillating pressure to the fiee stream 
dynamic pressure as high as 1. 4 per cent. It is feasible to present a plot oi pressure 
versus vehicle length for the boundary layer noise since the maximum dynamic pressure 
does not change appreciably when a liquid propellant booster is utilized. This is tiue 
because of the relationship between mass flow and engine thrust. Since the ratio oj. fiist 
stage thrust to total vehicle weight has crystallized around 1. 2 to 1. 5, and the flow rate- 
thrust relationship for liquid propellant boosters exists; there is only a small variation 
in the maximum dynamic pressure for any flight trajectory. Note that the distance in 
Figure 9 is measured from the leading edge of the vehicle. It was implicit y assumed 
that there would be no large projections or deviations from a smooth aerodynamic pro - 
file on the vehicle. Acoustic pressures can be magnified by a factor of 10 in the presence 
of an oscillating shock wave from an aerodynam ic discontinuity,. 

With only a knowledge of the approximate component location in respect to ve- 
hicle length, the curves may be utilized to determine each of the three pressures. Ob- 
viously, the most severe pressure should be utilized in the equation. Since equation 
(1) depends upon reference quantities these may be lumped as a constant. Further the 
component weight attenuation factor was defined to be (Appendix I) : 
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F = 


W 

w + w 

c 


where W = weight of structure 

W = component weight, 
c 

Substituting into equation (1) results in: 



where 



a constant depending upon the reference parameters 


W = panel weight for Type 1(a) , stiffener weight for Type 1(b) , and bulkhead 
weight for Type 1(c). 


In some cases a section of vehicle skin may not incorporate longitudinal stiffeners. 

When determining the panel weight, the radial distance should not exceed three times the 
distance between rings. This is because the dynamic characteristics of thin plates re- 
main essentially constant above aspect ratios of three. The acoustic pressures associ- 
ated with the reference vehicle environments are obtained from the predicted pressure 
at the appropriate station numbers. Although these pressures were measured during 
captive tests, it is anticipated that any inherent error associated with the acoustic pre- 
diction technique will appear in both the reference and new vehicle acoustic predictions 
with the same magnitude. Therefore, since ratios of pressure are involved, the error 
is considered to cancel. Thus, with adequate reference parameters, a knowledge of the 
component location with respect to vehicle length and the structural mass characteristics, 
the vibration response of the structure may be determined. If the reference environment 
represents the unloaded structure (i. e. , does not reflect the effect of component mass) 
then the environment obtained from the prediction techniques will correspond to the un- 
loaded structure of the new vehicle. This structural response will then represent an in- 
put to a component mounted on this structure. For this reason, the attenuation factor 
(F) must be incorporated into the equation to account for the component mass effect. In 
some instances during the preliminary development period, the weight of the basic struc- 
ture or component may not be known. In this case, the attenuation factor is considered 
as unity. This will result in an estimate of the component environment which may tend 
to be conservative. 
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Type II 


Now consider the Type II structure. This structure may be subdivided into 
two sections. These are: 

11(a) — Structural beams such as I beams, etc. The components mounted 
in this section would not primarily be effected by acoustic pressures but by rocket 
engine vibrations. 

H(b) — Rocket engine components. These components may again be sub- 
divided into three sections: 

11(b) -1: Combustion chamber section. This section includes the 
components mounted on the chamber dome or side case. 


11(b) -2: Turbo pump section. This section includes the components 
located on the propellant pumps. 

11(b) -3: Actuator assembly. The components located in this region 
are mounted on the actuator struts or actuator rods. 


The equation for predicting vibrations of these structural types is: 


G 

n 


G 



(NTV) 

(NTV) 


n 


r 


W 

r 


W 

n 


— v 
F 


(3) 


where G = the vibration response of the new structure. This is considered an input 
to the component mounted on this structure. 


G = the vibration response of a reference structure. This environment is 
determined by measured data. 


(NTV) = number of engines, thrust and exhaust velocity of rocket engine associated 
n with the stage under consideration in the new vehicle. 


(NTV) = number of engines, thrust and exhaust velocity of rocket engine associated 
with a reference vehicle. 

= the structural weight corresponding to G n 

= the structural weight corresponding to G^ 

= the weight of component to be mounted on 
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w 

n 


W + W 
n c 


= an attenuation factor which takes into account the effect of 
component mass loading 


Equation (3) applies to random rms composite values or to sinusoidal values. For the 
mean square spectral density value the total expression is simply squared. 


The geometry and mass characteristics of structural Type 11(a) are the most 
difficult to define; therefore, the problem must be approached in general terms. If the 
region under consideration is of beam and truss construction, the total structural weight 
(do not include the weight of components in the region) should be used to determine the 
parameter W n . Subsequent knowledge of separate component weights (W c ) would permit 
calculation of an individual attenuation factor (F) for each component in the region. This 
method is considered valid provided the reference vibration (G r ) is determined from an 
adequate number of measurements in this area and subsequent statistical analyses pro- 
vides confidence limits. With this criterion the vibration environment of any particular 
component located in this region may be specified with statistical certainty. 


W, 


The reference parameters 


may be lumped as a constant which simpli- 


(NTV) r 

fies the expression. Incorporating this constant and the expression for the attenuation 
factor into equation (3) results in: 


H(a) 




(NTV) 

n 


W + W 
n c 


(4) 


for 11(a) 

t 

for H(b) 


W = weight of entire beam area of new vehicle for 11(a) 
n 

W = weight of a single rocket engine section corresponding to the stage under 
n consideration for 11(b) 

For Type II (b) structure, clustering does not generally affect engine vibrations; there- 
fore, N r and N n are not considered. For a thorough understanding of equation (4) , a 
derivation is provided in Appendix n. Also, for type II structure, only one case is con- 
sidered the maximum. This case is the firing (either captive or inflight) of that 
particular stage under consideration. The vibration transmitted from other stage firings 
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C = 


C = 


1 


w 


(NTV) 


'\ 


W 


(TV). 




(i. e. , booster to upper stages) is negligible compared to that stage operation. Thus, 
with a knowledge of the necessary rocket engine parameters and structural mass charac- 
teristics, the vibratory environment of any particular component may be determined. 

An example problem dealing with this type of structure is provided in Appendix IV. 

Adequate reference parameters have been obtained from Saturn I Block I captive 
firing data. Currently, the Saturn I R&D Program has provided a total of 41 captive 
firings representing 2, 621 seconds of test time. These have resulted in 2, 708 vibration 
measurements representing 200 different locations throughout the booster. Approxi- 
mately 350 of these measurements have been utilized to define Type I structure and 725 
to define Type II structure. The remainder have been used to provide design modifica- 
tions for the operational Block II vehicles. Comprehensive statistical analyses were 
performed on these abundant data. With the results of this program, the respective 
environments of the eight structural types may be specified with statistical certainty. 

The statistical results, coupled with the necessary structural parameters, provide 
entirely satisfactory standards. The resultant environments and the associated param- 
eters corresponding to the unloaded vibratory response of the structural types is given 
in Figures 10 through 25. These environments are presented in both sinusoidal and 
random form and are to serve as reference criteria. The reference spectra are shaped 
according to the proportion of power within a specified bandwidth. The envelopes were 
constructed in a manner which will cover a slight shift of modal frequencies. Specifi- 
cally, large vehicles may be expected to possess a lower frequency content than smaller 
vehicles. 

The sinusoidal and random reference environments presented for each respective 
type of structure are not intended to be exact equivalencies. Since the present state of 
the art presents no valid random-sine correlations, the respective environments are 
considered satisfactory interpretations based upon the limited knowledge of random- 
sine relationships and successful design and testing results. 

Many sinusoidal testing techniques require both sweep and resonant dwell tests. 
The sinusoidal reference environments represent dwell criteria. For purposes of sine 
sweeping testing, the dwell level should be multiplied by an appropriate factor. 

With the accumulation of additional data, the standards may be slightly modified. 
The revision will, however provide a more exact definition of that particular structural 
type. Statistical methods of variance analysis are being applied to the Saturn I refer- 
ence structure to determine when a sufficient amount of data have been acquired. Type 
II structure has reached 95 per cent of the limit and Type I structure has reached 85 per 
cent of the limit. 
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CONCLUSIONS 


The methods and techniques presented herein will provide satisfactory estimates 
of the vibratory environment associated with any particular problem. These techniques 
are considered indicative of the present state-of-the-art and will permit satisfactory 
environmental estimates with only a few simple calculations. It is not imperative that 
the user thoroughly understand the philosophy behind the methods. He should, however, 
realize the necessary requirements and limitations. A summary of the principal limi- 
tations regarding these techniques is given below: 

1. The techniques do not apply to entire structure such as an entire rocket 
engine assembly or large structural members. They do, however, apply to components 
and other items of equipment mounted on these structures. 

2. These methods are not applicable for prediction of combined environments 
such as excitation due to propellant flow and engine vibration combined. In some remote 
cases the type of structure may not be clearly defined or excitation may be both acous- 
tically and mechanically induced. In these cases the two methods may be combined and 
a certain percentage assigned to each method depending upon degree of structural 
susceptibility to acoustic and mechanical excitation. Employment of the respective per- 
centages is left to the discretion of the user. Thus: 


where 


and 


The method is based upon the fact that similar structure has similar response character- 
istics (Appendix I page 21) * Precise frequency calculations tend to complicate the 
problem and it is doubtful that the results would justify the time and effort required. The 
calculation of the component natural frequency is, within itself, possibly a simple prob- 
lem. In the actual system, however, this frequency is coupled with the support structure 


VI 


X% G 2 +Y% G 2 


n 


II 


G = the resultant combined environment 
n t 

G = the environment obtained from using Type I methods 


G = the environment obtained from using Type n methods 
n II 


X% + Y% = 1. 0 or 100% 


3. This method does not provide frequency characteristics of the structure. 
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frequency characteristics. The determination of structure-component coupling effects 
is a highly complex requirement. Consequently, adequate frequency calculations and 
evaluation of these effects require considerable time and effort. 

A summary of the prediction equations is given below in outline form. The 
corresponding reference parameters are lumped as a constant in each of the respective 
equations. 


SUMMARY OF EQUATIONS FOR PREDICTING LOCALIZED 
VIBRATION ENVIRONMENTS OF A NEW VEHICLE 


t r ( to ) 


For Type I structure the constant represents ~ ; o (determined from Saturn I 

data) p r (lb/in > 


For Type n structure the constant represents^ ^ ^ ^ j/ S g C jj" (determined from 

Saturn I data) . 

For all structural types, is determined from Saturn I data. 


w imr 


Type I (a) Skin Panels 



where 

= panel weight 

G = Ref. Fig. 10 and 11 
r 

Type I (b) Ring Frames and Stringers. 



where 

W = stiffener weight 
n 

G = Ref. Fig. 12 and 13 
r 
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Type I (c) Bulkheads 


I (c)- 
G 

n^ 

G 

r 

I (c)- 
G 

n_ 

G 

r 

where 

W = 
n 

G = 
r 

Type II (a) 
G 

n_ 

G 

r 

where 

W 

n 

G - 
r 

Type n (b) 

n (b) 

G 

n_ 

G 

r 

where 

W 

n 


1 Forward Bulkheads 



2 Aft Bulkheads 



bulkhead weight 

Ref. Fig. 14, 15, 16, and 17. 


Beams 


= 9. 38 x 10 -4 


1 


(NTV) 


n 


(w n + w c ) 


= weight of entire beam area 
- Ref. Fig. 18 and 19 
Engine Components 
-1 Combustion Chamber Section 


= 7. 6 x 10 -4 


1 


(TV) 

n 

(W n + w c ) 


= combustion chamber + expansion nozzle weight 
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II (b) -2 Turbopump Section 


G 

n_ 

G 

r 

where 

W 

n 

n (b) 
G 

n_ 

G 

r 

where 

W 

n 

G 

r 


■\ 


= 6. 11 x 10" 4 


1 


(TV) 

n 

(W n + W c 


fuel + oxidizer pump weights 
3 Actuator Assembly 

(TV) 

n 

(W n + W c ) 


= 4. 87 x 10 


-4 


1 


weight of actuator rods and struts 
Ref. Fig. 20, 21, 22, 23, 24, and 25. 
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APPENDIX I 


DERIVATION OF EQUATION FOR PREDICTING THE VIBRATION 
ENVIRONMENT OF COMPONENTS MOUNTED ON SKIN STRUCTURE 


Any structure responding to oscillating forces possesses a vibrational power. 
For the case of skin structure the primary driving function is the impinging acoustic 
pressure which may also be defined in terms of power. Thus, a vibrational efficiency 
factor may be defined as: 


y = 


vib 


P 

ac 


( 1 ) 


This relationship stipulates that a certain amount of impinging acoustic power is trans- 
ferred into vibrational power. Equation (1) may be written for two separate vehicles as 



where the subscript r denotes a reference vehicle in which the vibratory environment 
has been defined by many measurements and n denotes a new vehicle of unknown environ- 
ments. For similar structure the above equations may be expressed as: 



Respectively, the average vibrational power and acoustic power may be expressed as: 
P vib = ^ ° 2 /cps W g Q (3) 

where G 2 / cps = power spectral density of the vibrating structure 

W = effective weight of the structure 
g^ = acceleration of gravity 
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and 


P = 


ac pc 
*a a 


(4) 


where p 2 = mean square acoustic pressure 

A = area over which the pressure is applied 

P a = mass density of the acoustic medium 

c = velocity of sound in the medium 
a 

Substituting equations (3) and (4) into equation (2) results in: 


(G 2 /cps) n W n (p a c a ) n P 2 A r Tn. 
(G 2 /cps) r W r (p a c a ) r P 2 A n = y r 


(5) 


Further, the transfer mechanism (y) is defined as: 


y = F 


2—i 


1 + 


m 


2Z 


a . 


where 

Z = the mechanical impedance per unit area of the structure (characteristic 
impedance) 

Z = the characteristic impedance of the acoustic driving function 

a 

For most practical cases Z « Z . Thus: 

1 am 


y 


a 


m 


Substituting equation ( 6) into equation (5) and rearranging the terms gives: 


( 6 ) 


P n (p a c a> 


«GVcp S ) n = (OVops) r ^f) r ^ p , 



(7) 


17 


The mechanical impedance of a structure may be expressed as: 


Z = 
m 


+ (Mw - — ) 
' co 


W 

S n 


C0‘ 


(V u n ) 


+ CO - 


n 


CO 


( 8 ) 


where D, M, and k are the effective damping, inertial, and elastic parameters of the 
system. Thus: 

(Z ) 
m n 


< Z m> 


w 

n 


% 

n 

W 

r 

r 2 f w n \ r\ 

(Ij“ n ) + (?- a, J 

% 

r 


Since it is assumed that similar structures possess similar dynamic characteristics 
(i. e. , similar damping and frequency content) : 


(Z ) W 

ra n. = __n 

(Z ) W 

m r r 


or 


< z m> 


n_ 


(Z ) 

m 


W 


n 


W 


( 9 ) 


Also the characteristic acoustic impedance may be defined as: 


Z = /3 p c 
a r K a a 


where 


(3 - a factor depending upon the geometry and coincident frequencies of a 
localized section of stiffened skin and is independent of the atmospheric 
impedance characteristics. 

Hence for similar structure: 


Substituting equations (9) and (10) into equation (7) gives: 


(G 2 / cps) ^ = (G 2 / cps) ^ 

W/A may be defined as weight per unit surface area. For isotropic panels .* 

W/A = p t (12) 

m 

where 



p = the weight density of the material 
m 

t = the thickness of die plate 

However, for corrugated or sandwich construction t is required to be an equivalent 
thickness. In regard to corrugation the weight over a particular surface (indicated by the 
sketch) is expressed as: 



. W - p t A = p t d S, 
m c m c s 


where 


(13) 


£ = surface length = Si 

s 
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Substituting equation (13) into equation (12) the expression for an equivalent "flat plate” 
thickness is shown to be: 


t = t 
eq c 



(14) 


This expression is only approximate since the corrugation exhibits a different set of 
dynamic properties than the isotropic plate. Also the geometry of the corrugation re- 
sults in variable elastic properties, within itself, at the boundaries ( anisotropicity) . 
However, since the resonant acceleration amplitude is primarily controlled by the sys- 
tem damping and mass properties, this technique is deemed sufficient. 


Sandwich plate construction may be handled in a similar manner. However, the 
damping properties of sandwich and homogeneous plates are significantly different and 
must be considered. From an apparent mass (acceleration impedance) analysis, the 
resonant response of thin plates is: 


F_ 

X 


pt An 

Sr, 


9 

4 


(15) 


where 

F 

= the apparent mass 

A 

rj = the structural loss factor 
For sandwich plates (refer to sketch for dimensions): 



h + h + h 
Pi P2 c 


where 

p = the density of the cover plates 
P 

P c = the core density 

Further, experimental investigations have indicated average damping values of: 
rj = . 044 (solid panels) 
tj = . 108 (sandwich panels) 


(16) 
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In order to establish equivalent plates, the resonant responses are required to be equal. 
Thus: 


Substituting the damping values and equation (16) into equation 17 results in: 


(17) 


t = 2. 5 
eq 


(p h + p h + p h ) 
Pl pi P2 P2 c C 

P m 


where 

p = the density of the equivalent homogeneous plate 
m 

If both the cover plates and core are of the same material then: 


t =2. 5 ( h + h + — h ) 
eq Pl P2 P m c 


(18) 


Further, if the cover plates are of the same thickness the above reduces to: 


t = 5 h + 2. 5 
eq p 


p h 
c c 


m 


(19) 


Generally, all rocket vehicles consist of aluminum skin which permits cancella- 
tion of the density terms. However, for the general case the densities will remain in 
the equation allowing calculation of the vibratory response for any type of skin material. 
Applying equation (12) to equation (11) results in: 


(G 2 /cps) n = (G 2 / cps) ^ 




( 20 ) 


If (G 2 /cps) r corresponds to an unloaded structure (i.e. , does not reflect the 
amplitude attenuation due to component mass) , then the environment determined by 
equation (20) would represent the unloaded response of the new structure. For this 
reason an attenuation factor is required to account for the component mass effect. Equa- 
tions (1) and (3) may be combined resulting in: 


7 p ac = 


1 

2ir 


G 2 /cps 


W 


§o 


( 21 ) 


This expression indicates that for a constant driving power the acceleration response 
decreases as the weight increases. Now consider an item of weight, W c , to be mounted 
on a structure in which the vibration environment has previously been defined. The 
effective force (acoustic pressures) driving this structure does not change. Hence, the 


vibration amplitudes will be decreased by a factor (F) of 


W 

W + W r 


or 


G 2 /cps (unloaded) = 


P vib 2 * f 
AfWg 0 
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and consequently 
G 2 /cps ( loaded) = 


P vib 2 * f 


W 


w + w. 


G 2 /cps (unloaded) 


(W + W c ) Afg 0 

Further the weight of the structure is an effective weight. This may be expressed as: 
W eff =IW act 

where I = a constant depending upon the system stiffness and mass coupling characteris- 
tics. This relation may be extended to include the added component weight by: 

W = I (W + W ) 
eff act c 

Substituting this equation into equation (3) does not alter the results since the constant 
cancels out 


Gycps (loaded) = 


4 


W 


(W + W c ) 


G 2 /cps (unloaded) 


( 22 ) 


A typical graph of this attenuation phenomena is shown in the sketch below: 



Incorporating equation (22) into equation (20) results in: 

W 


r P \ 2 f P 2 

(G 2 /cps) n = (G*/cps) r 1 n r r 


n 


W +W 
n c 


(23) 


For the composite or sinusoidal case the effective bandwidths cancel resulting in: 

V 


G = G l-B-lu! ^ 
11 r \pj VP. 


’ t /Aw 

n n J \ n 


W 


n 


+ W 


(24) 


The entire concept is based on the assumption that similar structure exhibits similar 
dynamic characteristics. Once the frequency and amplitude characteristics of a typical 
type of structure have been adequately defined, these results may be utilized as a basis 
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for predicting the vibration environment of like structure. The remaining factor to 
consider is the driving function. In the case discussed above the acoustic (or aerody- 
namic) pressures are the primary cause of vibration. These pressures must be pre- 
dicted for the new vehicle in order for equation (24) to be solved. 
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APPENDIX II 


DERIVATION OF EQUATION FOR PREDICTING 
VIBRATION ENVIRONMENT OF COMPONENTS MOUNTED 
ON STRUCTURE NOT SUSCEPTIBLE TO ACOUSTICS 


Dense structure such as heavy "I" beams are not affected by acoustic pressures. 
This type of structure is dynamically excited by rocket engine vibrations. Consider an 
engifte component subjected to a dynamic force that is random in nature. A certain 
amount of this force is "absorbed" by the component. The result of this absorption is 
exemplified as a dynamic response of the component which may be specified in terms of 
power. Hence, the mechanical efficiency factor is: 


y' 


vib 


mech 


where 


( 1 ) 


P = vibrational power 

vib 

P , = mechanical power 

mech 

This relation stipulates that a certain percentage of the potential mechanical power is 
"absorbed” as vibrational power. To further define the relationship, consider the ex- 
pressions for mechanical and vibrational power. 


p vib - h T 1 GVops w g o 

where 

W = the effective structural weight 

G 2 /cps = power spectral density of the vibrating structure 

g = acceleration due to gravity 
o 

Further: 

P = TV 
mech 

where 

T = thrust of the rocket engine 
V = exhaust velocity of the rocket engine 


( 2 ) 
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Substituting these expressions into equation (1) results in: 


Af W G 2 / cps g Q 
2tt f TV 


Assuming similar structure possess similar dynamic characteristics the mechanical 
efficiency factors may be considered equal. Thus: 


y' = y 
n 


f 

r 


(For similar structure) 

Where the subscripts denote new and reference vehicle, respectively or 


(G 2 / cps) ^ = (G 2 / cps) ^ 



This expression may be modified further by the mass attenuation factor (Appendix I) . 
Thus: 


(G/cps) n = (G 2 /cps) w 



(G 2 /cps) 



r\W +W 
' n c 


For the composite or sinusoidal case the effective bandwidths cancel, leaving: 


G 

n 




(T V ) 
n n 

(T V ) 
r r 


For predicting vibration environments of thrust frame areas composed of structural 
beams, the total mechanical power developed by that particular stage should be used. 

The reason for this is that the mass and geometry characteristics of thrust frame areas 
cannot be defined precisely enough to establish specific criteria. Therefore, the prob- 
lem must be approached in general terms. This dictates use of the total mechanical 
power and total thrust frame area weight. The reference vibration (G r ) should be based 
upon adequate statistical analyses of measurements taken at many locations throughout 
the thrust frame area. Thus, G r may be specified with statistical certainty and the 
general environmental prediction applied. The equation for this prediction is therefore: 
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(beams) 



(NTV) n 

(NTV) 


(3) 


where 


W = total weight of reference beam area 
r 

W = total weight of new beam area 
n 

W = weight of component to be located in this area 
c 

N = number of engines associated with that particular stage 

In the case of engine components, the number of engines is not necessary for 
determining the component environment. Investigations have indicated that clustering 
does not appreciably affect engine component vibrations. Thus: 


(engine components) 


G = G 
n r 



W 


w + w 

n c 


■) 


(TV) 

n 

(TV) 

r 


W 

r 


= basic weight of a single reference engine 


(4) 


W = basic weight of a single new engine 
n 

W = weight of component on the new engine 
c 

If the engine is separated into sections, then the basic weight is the weight of each indi- 
vidual section. These do not include component and hardware weights. It should be noted 
that the skeleton structure comprises approximately 80 per cent of the total engine 
weight. The requirements corresponding to equations (3) and (4) are that G r represent 
the unloaded structure (i. e. , does not reflect the effect of component mass loading) and 
that adequate measurements have been obtained to satisfactorily define the vibration 
environment. 
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APPENDIX IE 


EXAMPLE PROBLEM UTILIZING 
TYPE I METHODS 


Statement of Problem 

A component having a weight of 10 pounds is to be mounted on a 22 x 30 x . 15 
corrugated aluminum pjpiel as shown in the sketch below: 



The panel is located 800 inches above the engine nozzle exit plane of the second stage 
engines of the hypothetical vehicle configuration shown in enclosure (1). 

To obtain a reliable estimate of the components' vibrational environment, the 
following information concerning the overall vehicle configuration should be known: 


1. Total length of the vehicle configuration. 

2. Distance from engine nozzle exit plane of one stage to the engine nozzle 
exit plane of the preceeding stage for all stages of the configuration. 


3. Type and number of engines for the first and second stages. 


The type of structure for the given problem is in category I (a) skin structure, 
equation for predicting the environment of a component on this type of structure 
on page 12 . Thus: 



The- 
is found 
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where the material densities cancel because p r and p n are aluminum. 

Since the structure is susceptible to dynamic pressures, the maximum acoustic 
or aerodynamic pressure must be determined. Calculation of three pressures is neces- 
sary of which the highest pressure is used to determine the vibration environment. 

First, determine the overall acoustic pressure at the component location due to the 
captive firing of the second stage. From Figure 1, the second stage is equipped with 
three M-l engines. The equation for finding the resultant pressure produced by opera- 
tion of three M-l engines at a vehicle length of 800 inches is: 



From Figure 8, Pj is . 121 psi; therefore, the resultant pressure is: 
p^ = ^1T (. 128) = . 222 psi 

Next, determine the maximum overall acoustic pressure at the component loca- 
tion due to the firing of the booster stage. The distance of the component location from 
the booster engine nozzle exit plane is: 

1700 + 800 = 2500 inches 

From Figure 7, Pi = . 0435 psi 

Thus: 


P t = /^9 (.0435) = . 1305 psi 

Next, determine the maximum aerodynamic pressure. The distance of the 
component location to the leading edge of the vehicle is: 

1000 + 800 + 200 = 2000 inches 

The maximum aerodynamic pressure from Figure 9 is: 

p = . 1 psi 
ae 

From the selection of the most extreme pressure, 

p = . 222 psi 
n 

The equivalent thickness may be determined from Appendix I, equation (13). 
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Thus: 


(. 15) ( 3 + 3 + ^ -1 = e i 9 i inches 



Further, the panel weight is (from equation (12) , Appendix I): 
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W = p t dJi = (. 1) (. 15) (30) (3 + 4 + 3 + 4) 77 = 12. 6 lb 
m c s ll 


Substituting the calculated values into equation (1) results in: 


G 

n 

G 


r 


1. 22 





12. 6 


12 . 6+10 


-= 1. 06 


Now, the resultant environments may be obtained by multiplying the ratio 
reference environments. (Figs. 10 and 11). 


G 

n 


G 

r 


by the 


Hence, the vibratory environment in the given problem is: 



(Resonant Dwell 
Equivalent) 
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APPENDIX IV 


EXAMPLE PROBLEM UTILIZING 
TYPE II METHODS 


Statement of Problem 


Consider an item of equipment which weighs 30 lb. This component is to be 
mounted in the third stage frame area of the hypothetical vehicle shown in Enclosure (1). 
This area is composed of a network of structural beams. A compreshensive weight 
breakdown indicates a total thrust structure weight of 9, 000 lb. The problem is to 
determine the vibration environment of the component. 


The first step is to determine the type of structure where the component is to be 
located. From the definition of structural types indicated on page 8 , the problem ap- 
plies to type II (a) structure. The equation for determining the vibration environment 
of type H (a) is shown on page 13 . Thus: 


G 

n 

G 


r 



(NTV) n 

(NTV) r 


W 
r 

W 


n 


—i 

F 


9. 38 x 1(T 4 


(NTV) 


n 


(W +W ) 
' n c 


To solve this equation the following information is required: 


(1) the number of rocket engines associated with that particular stage 

(2) the thrust and exhaust velocity (ft/ sec) developed by the rocket engine. 


With this information: 


G 

= 9. 38 x 10~ 4 

Applying this factor to the sinusoidal and random reference environments shown in 
Figures 18 and 19, the resultant vibration environment of the component is: 

(3.2) (1.27) = 4. 06 G pk 
(.0028) (1.27) = .00356 In. D. A. Displ. 

(10.4) (1.27) = 13.2 G pk 


(6) (2 x 10 5 ) (1. 39 x 10 4 ) = 1 27 

(9,000 + 30) 
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(Resonant Dwell 
Equivalent) 



and for the random: 

(18.9) (1.27) = 24.0 Grms 

(. 18) ( 1. 27) 2 = .29 G 2 / cps 
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OBTAIN A KNOWLEDGE 
OF THE OVERALL 
VEHICLE CONFIGURATION 


VERIFICATION AND 
ANALYSIS OF THE 
DESIGNED COMPONENT 




DETERMINE THE TYPE AND 
WEIGHT (IF POSSIBLE) 

OF STRUCTURE WHERE 
COMPONENT IS TO BE 
LOCATED 


TEMPER WITH 
ENGINEERING 
JUDGEMENT 


«*■ 


TYPE I 

ESTIMATE CRITICAL 
PRESSURE (AERODYNAMIC 
OR ACOUSTIC) 


TYPE n 


DETERMINE THRUST AND 


EXHAUST VELOCITY OF 


ENGINES FOR THAT 


PARTICULAR STAGE 



SELECT 
CHOICE OF 
DESIGN 
METHOD 


USE KNOWN WEIGHT OF 


COMPONENT AND DETER- 


MINE THE VIBRATIONAL 


ENVIRONMENT 



ASSUME VIBRATIONAL 


ENVIRONMENT AND 


DETERMINE THE WEIGHT 


OF THE COMPONENT 



FIGURE 1 


FLOW CHART FOR DETERMINING THE VIBRATION ENVIRONMENT OF A 
NEW VEHICLE STRUCTURE 







MAXIMUM OVERALL RMS ACOUSTIC PRESSURE POUNDS/ (INCH) 



10 inches 

or less VEHICLE LENGTH (INCHES) 

(DISTANCE ABOVE ENGINE NOZZLE EXIT PLANE) 


FIGURE 2. ACOUSTIC PRESSURE VS VEHICLE LENGTH 165K H-l ENGINE 
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10,000 


10 Inches 
or less 


VEHICLE LENGTH (INCHES) 

(DISTANCE ABOVE ENGINE NOZZLE EXIT PLANE) 

FIGURE 3. ACOUSTIC PRESSURE VS VEHICLE LENGTH 188K H-l ENGINE 
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MAXIMUM OVERALL RMS ACOUSTIC PRESSURE POUNDS/(INCH) 



10 inches 

or less VEHICLE LENGTH (INCHES) 

(DISTANCE ABOVE ENGINE NOZZLE EXIT PLANE) 

FIGURE 4. ACOUSTIC PRESSURE VS VEHICLE LENGTH RL-10 ENGINE 
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MAXIMUM OVERALL RMS ACOUSTIC PRESSURE POUNDS/ (INCH) 



10 inches 

or less VEHICLE LENGTH (INCHES) 


(DISTANCE ABOVE ENGINE NOZZLE EXIT PLANE) 

FIGURE 5. ACOUSTIC PRESSURE VS VEHICLE LENGTH AGENA B&D ENGINE 
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MAXIMUM OVERALL RMS ACOUSTIC PRESSURE POUNDS/ (INCH) 



10 inches 

or less VEHICLE LENGTH (INCHES) 
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(DISTANCE ABOVE ENGINE NOZZLE EXIT PLANE) 

FIGURE 6. ACOUSTIC PRESSURE VS VEHICLE LENGTH J-2 ENGINE 


MAXIMUM OVERALL RMS ACOUSTIC PRESSURE POUNDS/ (INCH) 


10.0 


K0 



10 inches 
or less 


100C 


VEHICLE LENGTH (INCHES) 

(DISTANCE ABOVE ENGINE NOZZLE EXIT PLANE) 

FIGURE 7. ACOUSTIC PRESSURE VS VEHICLE LENGTH F-l ENGINE 


lO.OOff 
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MAXIMUM OVERALL RMS ACOUSTIC PRESSURE POUNDS/ (INCH) 



10 inches 
or less 


VEHICLE LENGTH (INCHES) 

(DISTANCE ABOVE ENGINE NOZZLE EXIT PLANE) 
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FIGURE 8. ACOUSTIC PRESSURE VS VEHICLE LENGTH M-l ENGINE 






RMS EFFECTIVE (5-2000C/S) PRESSURE (PSI) 



DISTANCE FROM LEADING EDGE OF VEHICLE 


(INCHES) 

FIGURE 9. BOUNDARY LAYER INDUCED NOISE VS VEHICLE LENGTH 


41 





00 (O o 




FIGURE 11. RANDOM REFERENCE ENVIRONMENT - TYPE 1(a) SKIN PANELS 
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FREQUENCY (CPS) 

FIGURE 12. SINUSOIDAL REFERENCE ENVIRONMENT - TYPE 1(a) SKIN STIFFENERS 



.001 


FREQUENCY (CPS) 


FIGURE 13. RANDOM REFERENCE ENVIRONMENT TYPE 1(b) SKIN STIFFENERS 
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FIGURE 14. SINUSOIDAL REFERENCE ENVIRONMENT TYjPE I(c)-i FORWARD BULKHEADS 



FIGURE 15. RANDOM REFERENCE ENVIRONMENT TYPE 1(c) -1 FORWARD BULKHEADS 
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.001 
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4 5 6 7 8 9 10 


1000 


FREQUENCY (CPS) 


FIGURE 17. RANDOM REFERENCE ENVIRONMENT TYPE 1(c) -2 AFT BULKHEADS 
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FIGURE 18. SINUSOIDAL REFERENCE ENVIRONMENT TYPE 11(a) BEAMS 



FIGURE 19. RANDOM REFERENCE ENVIRONMENT TYPE Ufa) BEAMS 
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FIGURE 20 



10 



FIGURE 21. RANDOM REFERENCE ENVIRONMENT TYPE 11(b) -1 COMBUSTION 
CHAMBER SECTION 
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FIGURE 23. RANDOM REFERENCE ENVIRONMENT TYPE U(b)-2 TURBOPUMP SECTION 


55 


CJl 





FIGURE 25. RANDOM REFERENCE ENVIRONMENT TYPE U(b)-3 ACTUATOR ASSEMBLY 
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